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Preface 


This  report  was  written  to  document  work  performed  in  the  area  of  verifica¬ 
tion  of  the  constitutive  model  used  in  the  nonlinear,  inaemental  structural  analy¬ 
sis  methodology.  The  work  was  sponsored  imder  funds  provided  to  the 
U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  by  the  Engineering 
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Computer-Aided  Analysis  of  Massive  Concrete  Structures  work  unit  of  the  Struc¬ 
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Director,  ITL.  Dr.  Robert  Hall,  Chief,  Structural  Analysis  Group,  SL,  was  the 
Program  Manager  for  the  Stmctural  Engineering  and  Conaete  R&D  Program. 
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The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or 
promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  for  the  use  of  such  commercial  products. 


Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  repOTt  can  be  converted  to  SI  imifs  as 
follows: 
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By 

To  Obtain 

degrees  Fahrenheit 

5/9 

degrees  Celsius  or  keMns^ 

feet 
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meters 
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25.4 

millimeters 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (mass) 

0.4535924 

kilograms 

1 '  To  obtain  Celsius  (C)  temperature  reading  from  Fahrenheit  (F)  reac 
To  obtain  Kelvin  (K)  readings,  use  K  =  (i 

lings,  use  the  following 

1  Introduction 


The  nonlinear,  inareroaital  structural  analysis  (NIS  A)  procedure  was  devel¬ 
oped  in  the  mid-1980s  concurrently  under  the  Cracking  of  Concrete  work  unit  of 
the  Rqtair,  Evaluation,  Maintenance,  and  Rdiabilitation  (REMR)  Program  and 
the  Melvin  Price  Locks  and  Dam  Project  The  NIS  A  procedure  was  developed  to 
update  a  methodology  used  since  the  late  1960s  to  evaluate  mass  conarete  struc¬ 
tures.  While  the  programs  used  for  evaluating  mass  conoete  structures  v^ch 
had  been  developed  by  Sandu,  William,  and  Raphael  (1967)  and  Wilson  (1968) 
had  been  used  successfully,  the  process  required  rqtdating  to  take  advantage  of 
the  inoreased  coroputational  capacity  of  computers  and  inqjroved  finitf.  ftlf.mftnt 
methods. 

The  efforts  resulted  in  the  development  of  a  procedure  to  be  used  for  the 
evaluation  of  massive  conarete  structures.  Following  its  use  on  the  Melvin  Price 
Project,  NIS  A  studies  were  performed  on  Lock  and  Dam  Nos.  4  and  5  on  the  Red 
River,  Olmsted  Locks  on  the  Ohio  Rivo",  and  the  McAlpine  Lock  Replacement 
Project,  also  on  the  Ohio  RivCT.  Formal  guidance  for  performing  a  NISA  was 
also  developed.  An  Engineer  Tedinical  Letter  (ETL  1 1 10-2-324),  “Special 
Design  Provisions  for  Massive  Concrete  Structures”  (U.S.  Army  Corps  of  Engi¬ 
neers  (USACE)  1990),  was  the  original  guidance  published  and  was  r^laced  by 
ETL  1110-2-365  (USACE  1994)  “Nonlinear,  Incremental  Stmctural  Analysis  of 
Massive  Conoete  Structures”  to  update  the  guidance  to  correspond  to  the  current 
state  of  the  art 

The  concurrent  effort  in  development  of  the  NISA  methodology  resulted  in 
many  procedural  issues  in  poforming  a  NISA  to  be  resolved  during  the  Melvin 
Price  Project  effort,  while  the  efforts  of  the  Cracking  of  Conoete  work  unit  con¬ 
centrated  on  the  development  of  a  constitutive  model  vriiidi  would  capture  the 
^e-dependent  characteristics  of  oeep,  shrinkage,  and  the  modulus  of  elasticity  of 
conoete.  In  addition,  a  smeared  oack  model  was  to  be  included  in  the  constitu¬ 
tive  model  to  allow  for  oack  prediction  in  structures.  The  genoal  purpose  finirf. 
element  package,  ABAQUS,  was  chosen  for  performing  NISA  studies  in  part 
because  uso  supplied  subroutines  could  be  used  for  constitutive  relationships 
After  a  decision  was  made  regarding  the  finite  element  code,  the  constitutive 
model  for  young  conoete  was  developed  by  ANATECH  Corporation  ntiliring 
their  smeared  oacking  and  oeep  model  used  for  evaluations  of  conoete  contain¬ 
ments  and  underground  waste  storage  tanks.  This  constitutive  model  is  contained 
in  the  proprietary  software  ANACAP-U  (ANATECH  1992b). 
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In  the  analysis  procedure  developed  by  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  and  ANATECH,  the  equations  in  the  constitutive 
model  were  calibrated  for  each  project  based  on  results  of  laboratory  tests  of  the 
actual  materials  expected  to  be  used  at  the  project  site.  A  review  of  the  equations 
within  the  constitutive  model  and  the  method  of  cahbrating  the  model  is  con¬ 
tained  in  WES  Technical  Report  ITL-95-8,  “Constitutive  Modeling  of  Concrete 
for  Massive  Concrete  Structures,  A  Simplified  Overview”  (Truman  and  Fehl 
1995).  This  report  details  the  parameters  required  to  model  elasticity,  creep,  and 
shrinkage  in  an  aging  model  and  how  to  modify  the  input  parameters  of  the 
constitutive  model  so  that  they  approximate  the  results  obtained  from  the  test 
data. 

Etespite  the  fact  that  the  NISA  process  was  established  based  on  past  experi¬ 
ences  in  performing  finite  element  analyses  of  massive  concrete  structures  and 
that  the  constitutive  model  was  being  calibrated  to  actual  test  results,  questions 
have  been  raised  concerning  the  validity  of  the  NISA  process.  Many  of  these 
questions  are  concerned  with  the  validity  of  using  laboratory  tests  to  capture 
behavior  that  is  occurring  within  a  massive  concrete  structure.  There  were  also 
concerns  about  the  validity  of  using  a  linear  viscoelastic  creep  model  based  on 
uniaxial  test  data  and  applying  it  to  all  states  of  stress. 

Due  to  the  concerns  regarding  the  NISA  procedure  and  the  constitutive  model, 
an  effort  was  undertaken  as  part  of  the  Structural  Engineering  and  Concrete  Pro¬ 
gram  work  unit  “Computer-Aided  Analysis  of  Massive  Concrete  Structures”  to 
validate  the  NISA  procedure.  This  effort  was  initiated  during  a  discussion 
between  a  group  of  engineers  and  scientists  with  extensive  experience  in  concrete 
testing  and  numerical  modeling  of  concrete.  The  discussions  centered  on  how  the 
NISA  methodology,  particularly  the  constitutive  model,  could  be  validated  and 
how  this  could  be  accomplished  using  existing  data  from  the  laboratory. 

The  result  of  the  discussions  led  to  a  proposed  plan  for  validating  the  constitu¬ 
tive  model  used  in  NISA.  The  following  validation  efforts  were  proposed: 

a.  Perform  analyses  of  slow  and  rapid  load  beam  tests  performed  for  the 
McAlpine  Lock  Replacement  Project  and  compare  analyses  results  with 
test  data. 

b.  Perform  creep,  shrinkage,  and  modulus  tests  using  mixture  proportions 
and  materials  from  the  McAlpine  project  to  check  the  capability  of 
repeating  the  results  of  the  original  set  of  tests. 

c.  Compare  model  predictions  of  creep  recovery  with  existing  test  data  to 
check  the  capability  of  the  numerical  model  to  produce  an  adequate 
representation  of  unloading  or  creep  recovery. 

d.  Perform  a  biaxial  creep  test  and  a  corresponding  finite  element  analysis  to 
determine  if  the  use  of  creep  curves  based  on  uniaxial  tests  provides 
accurate  predictions  for  other  states  of  stress. 
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e.  Test  a  representative  mass  concrete  placement  in  the  laboratory  with  con¬ 
trolled  boimdary  conditions.  Analysis  of  the  slab  should  be  performed 
prior  to  the  performance  of  the  test  to  establish  the  test  parameters. 
Multiple  lifts  are  necessary  and  cracking  should  occur. 

It  was  agreed  that  if  it  could  be  shown  that  the  constitutive  model  performed 
satisfactorily  in  the  above  tests,  then  the  model  could  be  used  with  a  high  level  of 
confidence.  It  was  also  recognized  that  the  results  of  the  proposed  plan  may 
show  that  adjustment  to  the  constitutive  model  is  necessary. 

This  report  addresses  items  a,  c,  and  d  of  the  above  plan.  Numerical  studies 
were  performed,  and  the  results  of  the  numerical  studies  are  compared  witii  the 
actual  data.  The  slow  load  and  rapid  load  beam  results  are  compared  with  results 
obtained  during  testing  perfomted  during  the  McAlpine  Lock  Replacement 
Project  The  creep  recovery  and  biaxial  numerical  studies  are  compared  against 
existing  test  data  rather  than  performing  a  new  test  This  varies  from  the  original 
proposal.  However,  since  a  literature  search  provided  reliable  test  data  on  biaxial, 
hydrostatic,  and  triaxial  states  of  stress,  the  study  was  expanded  to  include  these 
loading  conditions,  and  existing  data  were  used  for  the  comparison. 

Items  b  and  e  will  be  addressed  in  future  work. 
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2  Multiaxial  Creep 


Background 

Creep  is  the  time-dependent  deformation  of  concrete  subjected  to  a  sustained 
loading.  Creep  strain  is  the  total  strain  of  a  loaded  specimen  minus  the  elastic 
strain  that  occurs  upon  loading  and  the  shrinkage  and  thermally  induced  strains 
determined  from  control  specimens. 

The  majority  of  creep  tests  are  uniaxial  compressive  tests  conducted  in  accor¬ 
dance  with  American  Society  for  Testing  and  Materials  (ASTM)  Standard  C 
512-87  (1992a),  which  requires  that  the  maximum  test  load  be  no  more  than 
40  percent  of  the  compressive  strength  of  the  concrete  at  the  age  of  loading.  At 
levels  of  stress  below  this  value,  creep  is  assumed  to  be  a  linear  function  of 
stress.  At  higher  levels  of  stress,  creep  can  be  highly  nonlinear.  The  ANACAP-U 
model  uses  this  assumption  of  creep  linearity  to  predict  creep  strain  under  all 
states  of  stress.  Since  maximum  stress  in  mass  concrete  is  limited  by  the  tensile 
capacity  of  the  concrete  and  is  normally  significantly  less  than  40  percent  of  the 
compressive  strength,  this  assumption  of  linearity  is  applicable  for  uniaxial 
loadings,  both  compressive  and  tensile.  However,  there  is  much  less  data  avail¬ 
able  on  ftinhiaxial  creep  testing,  and  available  test  data  are  often  inconsistent 

Multiaxial  creep  tests  cited  in  American  Concrete  Institute  (ACI)  publica¬ 
tion  SP-27  (1970)  resulted  in  conflicting  conclusions  about  the  effects  of  state  of 
stress  on  creep.  However,  many  of  the  differences  were  attributed  to  differences 
in  environmental  conditions,  with  drying  being  the  most  influential  factor.  This 
is  a  logical  conclusion  since  creep  is  commonly  aknowleged  to  be  heavily 
dependent  on  moisture  migration.  The  ANACAP  aging  viscoelastic  model  was 
designed  for  use  on  mass  concrete,  in  which  drying  only  occurs  near  the  surface 
and  the  majority  of  the  concrete  is  isolated  from  moisture  loss.  To  accurately 
model  the  moisture  conditions  of  mass  concrete,  only  tests  on  sealed  specimens 
were  used  to  verify  ANACAP  multiaxial  predictions. 
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Description  of  Verification  Tests 


Creep  test  results  used  in  the  multiaxial  aeep  verification  were  reported  by 
York,  Kennedy,  and  Perry  (1970)  and  McDonald  (1973).  In  these  tests,  uniaxial, 
biaxial,  and  triaxial  aeep  tests  were  conducted  on  6-in.^  by  18-in.  or  6-in.  by 
16-in.  cylinders,  and  unloaded  con:q}anion  cylindCTS  were  monitored  for  shrinkage 
strains.  Test  cylinders  were  cured  using  two  different  processes;  as-cast  or  air 
dried.  As-cast  cylinders  were  kept  in  a  fog  room  for  1  to  2  days,  then  sealed 
against  moisture  loss  with  two  coats  of  epoxy  and  enclosed  in  copper  sleeves. 
Two  temperatures  at  loading  were  also  considered  in  each  series  of  tests:  75  °F 
and  150  °F  in  the  tests  by  York,  Kennedy,  and  Pery  (1970),  73  °F  and  150  °F  in 
the  tests  by  McDonald  (1973).  Cylinders  were  loaded  at  90  days  after  casting, 
and  loadings  were  maintained  for  1  year.  Since  conditions  in  the  interior  of  a 
mass  concrete  structure  could  be  expected  to  be  similar  to  those  of  the  as-cast  73 
or  75  °F  cylinders  at  90  days  after  placement,  only  the  tests  that  wcare  performed 
under  as-cast  enuring  conditions  were  modeled  in  the  finite  element  analyses. 

In  the  tests,  axial  loads  were  applied  by  a  hydraulic  ram  and  monitored  by  a 
load  cell.  Radial  loads  were  ^plied  by  hydraulic  pressure  through  oil  contained 
by  a  steel  pressure  jaclrat.  A  control  system  automatically  regulated  pressures  to 
within  5  percent  of  the  design  pressures  excq)t  in  the  event  of  a  power  failure  or  a 
major  oil  leak. 

Conq)ressive  strength,  elastic  modulus  E,  and  Poisson’s  ratio  at  loading,  and 
design  loadings  for  the  as-cast  test  C5dinders  are  listed  in  Table  1.  All  tests 
designated  by  a  single  letter  and  a  number  (for  exanq>le,  B-7)  are  by  York, 
Kennedy,  and  Pory  (1970). 


Description  of  FE  Model 

The  tests  were  simulated  using  the  finite  element  (FE)  program  ABAQUS  and 
the  aging  viscoelastic  material  model  ANACAP-U  developed  by  ANATECH 
Corporation.  A  one-quartCT  symmetric  section  of  a  cylinder  was  modeled  using 
ei^t-node  axisymmetric  elem^ts,  and  the  appropriate  uniform  pressures  were 
^plied  at  the  top  and  outer  side  of  the  grid.  The  grid  used  in  the  FE  analyses  is 
shown  in  Figure  1. 

Test  E-39,  a  600-psi  imiaxial  test,  was  used  to  calibrate  the  creep  curve  in 
ANACAP-U.  Due  to  lack  of  data,  no  atten:q)t  was  made  to  model  the  aging 
modulus.  Although  dianges  in  modulus  are  used  to  translate  the  creep  curve  in 
time  (ANATECH  1992a),  the  increase  in  modulus  after  90  days  is  relatively 
small.  Since  these  test  simulations  were  begun  at  an  initial  age  of  90  days,  errors 
due  to  using  a  modulus  curve  based  on  data  from  other  conaetes  should  be 


^  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  of  measurement  is  found 
on  page  vi 
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Table  1 

Test  and  Cylinder  Parametei 

90>Day 
Compressive 
Strength,  psi 

Elastic  Parameters 

Design  Pressure,  psi 

Test 

E,psix1(l' 

Poisson’s 

Ratio 

Axial 

Radial 

B-7 

6,110 

5.66 

0.243 

2.400 

0 

B41 

6,110 

7.14 

0.239 

1,200 

2.400 

016 

6,430 

2.04 

0.409 

1,200 

1,200 

023 

6,430 

6.54 

0.244 

2,400 

600 

026 

6.500 

5.91 

0.273 

3.600 

1,200 

031 

6,500 

3.11 

0.358 

3,600 

3,600 

E-5 

7,290 

2.78 

0.381 

600 

600 

E-39 

7,290 

6.04 

0.293 

600 

o 

F-9 

7.410 

12,55 

0.017 

2,400 

2,400 

F-13 

7,410 

6.13 

0.264 

0 

600 

035 

6,460 

4.91 

0.248 

600 

3,600 

H-22 

6,330 

4.51 

0.209 

0 

3,600 

McDonald  Axial 

7,660 

5.50 

2.400 

0 

McDonald 

Triaxial 

7,660 

5.50 

1 

2,400 

600 

small.  The  fit  to  the  test  E-39  strains  was  based  upon  the  McAlpine  exterior 
mixture  creep  and  modulus  curves  (Bombich,  Neeley,  and  Garner.)^  The  creep 
curve  was  then  modified  by  means  of  the  creep  factor  until  predicted  aeep 
strains  closely  matched  test  strains.  Predicted  and  test  axial  strains  are  shown  in 
Figure  2.  The  2,400-psi  axial  test,  B-7,  was  modeled  to  verify  this  fit  For  this 
analysis,  predicted  axial  creep  strains  wwe  94.4  percent  of  test  strains  on  day 
360.  Predicted  and  test  axial  creep  strains  are  shown  in  Figure  3.  These  axial  test 
gimnlatinns  produced  satisfactory  results,  and  the  calibration  was  then  used  in  the 
simulation  of  all  tests  by  York,  Kennedy,  and  Perry  (1970). 
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^  Bombich,  A.  A.,  Neeley,  B.  D.,  and  Gamer,  S,  B.  “Concrete  mixture  selection  and  characteriza¬ 
tion  study  McAlpine  locks  replacement,  Ohio  River,”  (technical  report  in  preparation),  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


Chapter  2  Muttiaxial  Creep 


Creep  Strain  (micro  in/in) 


Rgure  1.  FE  grid 


Figure  2.  Test  and  predicted  strains,  E-39 


Chapter  2  MuKlaxial  Creep 


Rgure  3.  Test  and  predicted  strains,  B-7 


FE  Analyses  Results 

Only  two  biaxial  tests  are  listed  in  Table  1:  test  F-13,  with  a  600-psi  radial 
pressure  and  test  H-22,  with  a  3,600-psi  radial  pressure.  Test  creep  strains  and 
creep  strains  predicted  by  the  FE  analysis  for  the  600  psi  test  compare  well  in 
both  the  radial  and  axial  directions,  and  are  shown  in  Figure  4.  Predictions  for 
the  3,600-psi  test  could  be  expected  to  be  much  different  from  test  strains  for  two 
reasons.  First,  the  Von  Mises  stress  of  3,600  psi  is  approaching  a  limiting  value 
at  90  days  and  could  result  in  a  shear  failure  if  predicted  compressive  strength 
was  significantly  less  than  test  strength.  Secondly,  the  creep  model  in  the 
ANACAP-U  is  designed  to  predict  creep  for  loadings  of  up  to  0.4/'^,  or  approxi¬ 
mately  2,500  psi  at  day  90.  In  this  range,  creep  is  assumed  to  be  linear,  and  a 
single  creep  curve  can  be  used  to  define  creep  under  any  loading  and  translated  in 
tinte  to  define  creep  at  any  age.  At  higher  loadings  this  linear  creep  relationship  is 
not  ^plicable.  To  ensure  that  failure  did  not  occur,  the  3-day  compressive 
strength  was  increased  from  1,260  psi,  the  Me  Alpine  value,  to  4,260  psi  for  the 
analysis  modeling  this  test  Analysis  predictions  for  this  test  were  much  lower 
than  test  creep  strains  (Figure  5)  due  to  the  nonlinearity  of  creep  at  this  stress 
level.  The  assumption  of  the  linearity  of  creep  is  valid  for  mass  concretes,  where 
high  compressive  stresses  are  unusual,  and  the  stress  range  of  interest  is  normally 
limited  by  the  tensile  capacity  of  the  concrete,  but  may  result  in  errors  in  pre¬ 
dicted  creep  strains  if  the  ANACAP-U  aging  viscoelastic  model  used  in  these 
analyses  is  used  to  model  structural  concrete,  which  can  experience  high 
compressive  stresses. 
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Rgure  4.  Test  and  predicted  strains,  600-psi  biaxial  loading  (F-1 3) 


Rgure  5.  Test  and  predicted  strains,  3,600-psi  biaxial  load 


Chapter  2  Muitiaxiai  Creep 


Predicted  creep  strains  for  the  triaxial  test  simulations  agreed  well  with  test 
strains  in  tests  C-23  and  D-26,  in  which  radial  stresses  were  low  when  compared 
with  axial  stresses.  Predicted  axial  creep  strain  for  the  C-23  test  on  day  360  was 
only  13  percent  greater  than  test  axial  creep  strain  at  that  time.  Predicted  axial 
and  radial  creep  strains  for  test  specimen  D-26  were  6  and  13  percent  greater  than 
the  respective  test  values  on  day  360.  Predicted  and  test  strains  for  these  two 
tests  are  shown  in  Figures  6  and  7. 
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Rgure  6.  Test  and  predicted  strains,  2,400-psi  axial  pressure,  600-psi  radial 
pressure  (C-23) 


For  load  cases  in  which  stresses  were  hydrostatic,  significant  differences  were 
observed  between  the  test  and  predicted  strains.  For  all  hydrostatic  loadings,  test 
radial  creep  strains  exceeded  test  axial  creep  strains.  This  was  in  part  explained 
by  firiction  losses  in  the  axial  loading  mechanism  that  resulted  in  actual  pressures 
3  to  10  percent  below  design  values.  However,  radial  strains  exceeded  axial 
strains  by  more  than  10  percent  in  all  hydrostatic  tests.  In  test  E-5,  with  a  600-psi 
hydrostatic  loading  (Figure  8),  predicted  creep  strains  were  10  percent  lower  than 
test  radial  creep  strains  and  9  percent  higher  than  text  axial  creep  strains  at 
120  days  after  loading.  This  difference  is  partly  attributable  to  the  fact  that  the 
design  loads  used  in  the  numerical  test  simulations  were  slightly  different  from 
the  actual  loadings  due  to  the  reduced  axial  loading.  However,  by  day  360  the 
predicted  creep  strains  were  23.2  percent  greater  than  the  radial  test  creep  strain 
and  90.2  percent  greater  than  the  axial  test  creep  strain,  which  decreased  after 
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day  120.  The  fact  that  the  test  gauges  indicated  zero  or  decreasing  creep  between 
120  and  360  days  would  seem  to  indicate  a  loss  of  load,  but  none  was  reported. 
In  the  2,400-psi  hydrostatic  test  (Figure  9),  predicted  axial  creep  strains  were 
close  to  the  test  strains,  but  predicted  radial  creep  strains  were  only  70  to  76  per¬ 
cent  of  the  test  strains.  For  the  3,600-psi  hydrostatic  test  (Figure  10),  predicted 
radial  creep  strains  were  only  about  75  to  79  percent  of  test  strains,  while  pre¬ 
dicted  axial  creep  strains  were  82  to  90  percent  of  test  strains.  Predictions  for 
this  test  are  better  than  expected  since  the  loading  is  outside  of  the  range  for 
linear  creep.  The  level  of  error  seen  in  these  predictions,  10  to  30  percent,  may 
be  an  acceptable  level  of  error  for  creep,  which  is  only  one  of  the  many  factors 
that  affect  stresses  and  strains  in  mass  concrete.  In  the  1,200-psi  hydrostatic  test 
(Figure  1 1),  radial  pressure  was  lost  after  90  min.  Predicted  axial  creep  strain  for 
this  test  was  47.8  percent  greater  than  the  test  strain  on  day  360.  In  a  separate 
analysis,  a  l,2(K)-psi  uniaxial  test  was  simulated  although  no  uniaxial  test  data 
was  available  for  this  loading.  The  resulting  axial  creep  strains  were  exactly 
equal  to  the  strains  predicted  for  the  C-16  loading.  In  other  words,  applying  the 
radial  stress  for  a  period  of  0.0625  days  (90  min)  had  no  effect  on  the  predicted 
axial  strains. 


Rgure  9.  Test  and  predicted  strains,  2,400-psi  hydrostatic  loading  (F-9) 

The  largest  differences  in  predicted  and  test  strains  occurred  in  tests  B-41 
(Figure  12)  and  G-35  (Figure  13).  In  test  B-41,  with  a  2,400-psi  radial  pressure 
and  1,200-psi  axial  pressure,  test  axial  creep  strains  were  much  larger  than  test 
radial  creep  strains.  As  would  be  expected,  ANACAP-U  predicted  larger  radial 
creep  strains  and  very  little  axial  creep  strain.  No  explanation  was  provided  for 
this  apparent  reversal  of  test  creep  strains,  although  a  radial  pressure  drop  was 
reported  on  day  21.  Reported  total  radial  strains  dropped  between  day  0.25  and 
day  1,  on  day  6,  and  on  day  21,  possibly  indicating  drops  in  radial  pressure. 
Since  pressure  histories  were  not  reported,  accurate  simulation  of  pressure  drops 
in  the  FE  analyses  was  not  possible.  Specimen  G-35  was  loaded  with  a  600-psi 
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Rgure  11.  Test  and  predicted  creep  strains  for  1 ,200-psi  hydrostatic  loading 
vwth  loss  of  radial  load  after  90  min  {C-16) 
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Creep  Strain  (micro  in/in) 


Rgure  12.  Test  and  predicted  strains  for  1 ,200-psi  axial  loading,  2,400-psi 
radial  loading  (B-41) 


Figure  1 3.  Test  and  predicted  creep  strsuns  for  600-psi  axial  loading,  3,600-psi 
radial  loading  (G-35) 
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axial  stress  and  3,600-psi  radial  stress.  Although  the  Von  Mises  stress  for  this 
condition  was  only  3,000  psi,  the  analysis  3-day  compressive  strength  was 
increased  to  4,260  psi  to  ensure  that  shear  failure  did  not  occur.  Test  radial  creep 
strains  for  this  specimen  were  75.8  percent  greater  than  ANACAP-U  radial  creep 
strains  on  day  360,  and  test  axial  strains  were  65  percent  greater  than 
ANACAP-U  strains.  However,  3,600-psi  is  approximately  56  percent  of  the 
90-day  compressive  strength  for  specimen  G-35  and  is  outside  the  range  of 
stresses  for  which  creep  can  be  considered  to  be  linear 

Two  tests  by  McDonald  (1973)  were  also  modeled.  The  2,400-psi  axial  test 
was  used  to  calibrate  the  material  model,  and  a  triaxial  test  with  2,400-psi  axial 
pressure  and  600-psi  radial  pressure  was  analyzed  using  the  calibrated  model. 
Predicted  axial  strain  for  the  triaxial  FE  analysis  was  13  percent  greater  than  test 
axial  strain  on  day  360.  Analysis  radial  strains  were  slightly  higher  than  test 
radial  strains,  which  were  close  to  zero.  Calculated  and  test  sttains  for  the  two 
tests  are  compared  in  Figures  14  and  15. 


Rgure  14.  Test  and  predicted  strains,  2,400-psi  axial  creep  test,  McDonald 
main  modulus  concrete 


Conclusions 


The  creep  model  in  the  ANACAP-U  viscoelastic  aging  concrete  model  is 
based  on  the  linearity  of  creep  for  loadings  in  the  range  of  0  to  40  percent  of  the 
compressive  strength  of  the  concrete.  This  assumption  is  certainly  applicable  to 
the  prediction  of  stress  and  strain  in  mass  concrete  structures,  which  normally 
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Rgure  15.  Test  and  predicted  creep  strains,  triaxial  test,  McDonald  main 
modulus  concrete 


experience  loadings  much  lower  than  the  40  percent  upper  limit.  Based  on  the 
multiaxial  test  data  of  McDonald  (1973)  and  York,  Kennedy,  and  Perry  (1970), 
calibration  of  the  model  to  uniaxial  data  yields  acceptable  predictions  for  biaxial 
and  triaxial  states  of  stress  that  fall  within  the  range  for  which  the  linear  creep 
model  was  intended.  At  loadings  greater  than  40  percent  of  the  concrete 
compressive  strength,  creep,  as  expected,  was  highly  nonlinear  and  strains 
predicted  by  the  model  were  much  lower  than  test  strains. 
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3  Creep  Strain  Recovery 


Background 

As  discussed,  the  aeq)  model  in  the  ANACAP-U  viscoelastic  aging  model 
relies  on  the  linearity  of  acep  undo-  all  states  of  stress.  This  means  that  the  samp 
creep  curve  used  for  loading  is  also  used  to  predict  change  in  cxesp  strain  due  to 
unloading.  Howeva:,  it  has  been  observed  in  WES  aeep  tests  in  \nirich  the 
loading  was  lost  due  to  leaks  in  the  hydraulic  system  or  power  failures,  that  the 
amount  of  strain  recova-ed  iqton  imloading  seemed  mudi  smalla:  than  that  vshich 
might  be  expected  for  such  tests.  A  7-day  creq)  test  in  \shich  the  loading  was  lost 
on  day  8  is  shown  in  Figure  16.  This  test  was  not  modeled  in  this  SCTies  of 
verification  tests  due  to  lack  of  necessary  test  data,  but  the  small  amnnnt  of  areqp 
strain  recovery  can  be  observed  in  the  figure.  Tha  may  be  due  in  part  to  the  fact 
that  the  conaete  was  a  high  fly-ash  mixture  that  was  not  expected  to  mature  until 
q>proximately  90  days,  so  any  loading  occurring  at  early  ages  could  result  in  a 
permanait  set  that  could  not  be  predicted  by  the  ANACAP-U  model. 

The  aeq)  tests  described  in  the  previous  chapto:  (York,  Kramedy,  and  Pmy 
1970)  woe  conducted  on  mature  conarete  and  ware  monitored  for  90  days  after 
the  loading  was  removed.  These  tests  have  bear  used  to  test  the  accuracy  of  the 
ANACAP-U  creep  strain  recovery  predictions.  In  each  of  the  tests,  the  aeep 
cylindo:  was  loaded  for  1  year,  at  whidi  time  the  load  was  removed,  and  strains 
were  monitored  for  an  additional  3  months  afto'  unloading  to  provide  information 
on  areq)  strain  recovoy. 


FE  Modeling  of  Creep  and  Creep  Recovery 

For  stresses  less  than  40  percent  of  concrete  stroigth,  concrete  is  commonly 
treated  as  an  aging  viscoelastic  material.  Creep  in  this  stress  range  is  assiuned  to 
be  linear.  This  means  that  at  any  point  in  time  strain  is  proportional  to  stress,  and 
CTeep  is  conpletely  defined  by  the  creep  compliance  curve.  The  oreep  compli¬ 
ance  curve,  as  defined  for  a  mature  cono-ete  by  ASTM  C  5 12-87  (ASTM  1992a), 
is  the  strain  per  unit  of  stress  resulting  after  all  shrinkage  and  thomal  strains  and 
the  initial  elastic  strain  have  been  subtracted  fix»m  the  total  strain  of  a  loaded 
specimen.  For  use  in  the  aging  viscoelastic  model,  in  whidi  modulus  and  there¬ 
fore  elastic  strain  dianges  with  time,  this  definition  of  spedfic  aeq)  has  beoi 
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Rgure  16.  7-day  creep  test  with  loading  lost  on  day  8  (Gamer,  Hammons,  and 
Bombich  1991) 

modified  by  subtracting  the  changing  elastic  strain  due  to  a  given  loading  from 
the  total  strain. 

The  assumption  of  linearity  allows  the  application  of  the  principle  of 
superposition  (Bazant  1988).  Using  this  principle,  current  creep  strain  can  be 
obtained  as  a  hereditary  integral  over  the  stress  history,  and  creep  is  recoverable 
to  the  extent  that  creep  recovery  is  limited  by  aging.  A  typical  creep  compliance 
curve  is  shown  in  Figure  17  (Truman  and  Fehl  1995).  In  this  figure,  creep  com¬ 
pliance  for  a  non-aging  material  is  given  by  /(/-to),  and  creep  con^)liance  for  an 
aging  material  is  given  by  C(/-to).  A  discussion  of  the  extension  of  the  uniaxial 
creep  relationship  to  the  multiaxial  creep  equations  used  in  ANACAP-U  may  be 
found  in  the  ANACAP-U  theory  manual  (ANATECH  Research  Corporation 
1992a). 

One  of  the  limitations  of  the  principle  of  superposition  discussed  by  Bazant 
(1988)  is  that  it  does  not  provide  accurate  predictions  for  a  decrease  in  the  strain 
magnitude  due  to  unloading.  However,  in  practice  superposition  is  commonly 
used  to  predict  strains  during  unloading  even  though  this  may  result  in  some  loss 
of  accuracy.  Also,  if  stresses  are  changed  very  slowly,  the  use  of  the  principle  of 
linear  superposition  is  generally  considered  to  predict  creep  strains  due  to  stress 
fluctuations  satisfactorily  (Bazant  1988). 

In  the  FE  analyses,  the  effects  of  creep  strain  recovery  were  simulated  by 
extending  the  analyses  (discussed  in  the  previous  chapter)  for  90  days  after 
removal  of  the  loadings.  For  purposes  of  comparison  between  predicted  and  test 
curves,  creep  strains  were  set  equal  to  zero  immediately  after  unloading  in  all 
curves.  Only  changes  in  strain  after  unloading  were  considered. 
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Rgure  17.  Typical  specific  creep  curves,  aging  and  non-aging  materiats 
(where  E=  modulus  of  elasticity) 


Analyses  Results 

Comparisons  between  predicted  strains  after  unloading  and  test  curves  for  the 
imiaxial  and  biaxial  tests  are  shown  in  Figure  18,  and  those  for  the  hydrostatic 
tests  are  shown  in  Figure  19.  In  the  uniaxial  and  biaxial  plots,  predicted  creep 
recovery  strains  in  the  direction  of  loading  were  between  200  and  300  percent  of 
test  strains.  For  the  hydrostatic  tests,  predicted  creep  strain  recovery  was  close  to 
three  times  the  test  creep  strain  recovery  in  each  of  the  three  tests  modeled. 

Predicted  and  test  creep  recovery  strains  for  three  of  the  triaxial  tests  are 
shown  in  Figure  20.  In  tests  C-23  and  D-26,  radial  pressures  (6(X)  and 
1,200  psi  respectively),  were  small  when  compared  with  axial  pressures  (2,400 
and  3,600  psi,  respectively)  and  both  test  and  predicted  strain  recoveries  were 
small.  Predicted  axial  creep  recovery  for  these  cylinders  was  greater  than 
200  percent  of  the  test  values.  In  test  G-35,  axial  pressure  (600  psi)  was  small 
when  compared  with  radial  pressure  (3,600  psi),  and  predicted  radial  creep 
recovery  was  greater  than  200  percent  of  the  test  values. 

In  Figure  21,  predicted  and  test  axial  creep  recovery  strains  for  test  E-39,  the 
600-psi  umaxial  test,  are  plotted  on  a  positive  scale  with  test  and  predicted  creep 
strains  under  loading  for  purposes  of  conq)arison.  Predicted  and  test  creep  strains 
at  80  days  after  loading  and  predicted  and  test  creep  recovery  strains  at  80  days 
after  unloading  are  given  in  Table  2.  Predicted  axial  creep  recovery  strains  were 
similar  to  predicted  axial  strains  under  loading,  as  would  be  expected  for  a  mature 
concrete  and  a  linear  creep  model,  while  test  axial  creep  recovery  strains  were 
only  about  one-third  of  the  strains  predicted  under  loading. 
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Figure  1 9.  Hydrostatic  test  and  predicted  creep  recovery 


Figure  20.  Test  and  predicted  creep  recovery,  triaxial  tests 


Rgure  21 .  Test  and  predicted  creep  and  creep  recovery  strains,  E-39 


This  plot  demonstrates  the  shortcomings  of  the  linear  creep  model  in  cases  of 
dramatic  unloading. 


Conclusions 

The  large  differences  between  test  and  predicted  creep  recovery  demonstrated 
in  the  figures  above  occur  under  circumstances  for  which  the  principle  of  super¬ 
position  could  not  be  expected  to  provide  accurate  predictions.  However,  these 
conditions  do  not  typically  occur  in  a  thermal  stress  analysis  of  young  concrete. 
Strains  in  mass  concrete  are  normally  induced  by  relatively  slowly  changing 
loadings.  Thermal  loadings  are  usually  provided  by  the  average  ambient  tempera¬ 
ture  curve  and  by  the  adiabatic  temperature  rise  of  the  material.  For  a  thick  wall 
placed  during  the  warmest  part  of  the  year,  tenqieratures  near  the  center  rise 
significantly  above  ambient  air  ten:q)eratures  due  to  the  heat  generated  by  the 
hydration  of  the  cementious  materials,  while  siuface  temperatures  remain  much 
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closer  to  the  ambient  air  temperature.  This  process  creates  a  temperature  gradient 
between  the  outer  surface  and  the  center  of  the  wall,  resulting  in  expansion  at 
the  center  and  contraction  at  the  surface  as  ambient  temperatures  fall.  These 
surface  contractions  are  restrained  by  the  large  mass  of  heated  interior  concrete 
and  by  the  thermal  gradient,  resulting  in  tensile  stresses  at  the  surface.  These 
surface  tensile  stresses  are  somewhat  relieved  by  the  reciprocal  effect  of  creep, 
stress  relaxation.  Under  the  conditions  described,  maximum  tensile  stresses 
should  occur  near  the  minimum  of  the  ambient  temperature  curve.  Surface 
stresses  will  then  slowly  decrease  as  the  average  ambient  temperature  rises.  The 
ANACAP-U  model  should  be  able  to  represent  creep  and  stress  relaxation  fairly 
accurately  under  the  conditions  described.  However,  creep  predictions  may 
become  less  accurate  for  a  sudden  drop  in  loading  since  creep  recovery  will  be 
overpredicted. 

At  present  there  is  no  consistent  data  to  justify  using  anything  other  than  a 
viscoelastic  creep  model  in  cases  of  unloading.  At  least  one  series  of  creep  tests 
in  which  the  loading  is  gradually  reduced  with  time  and  one  series  in  which  the 
loading  is  reduced  in  a  series  of  sudden  steps  would  provide  data  for  determining 
the  level  of  accuracy  of  the  current  model  under  these  conditions,  and  possibly 
provide  data  to  initiate  the  consideration  of  modifications  to  the  method  of 
predicting  creep  recovery. 
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4  Slow  Load  Beam  Test 


Introduction 

The  slow  load  beam  test  is  used  to  determine  the  ultimate  failure  strain  of  the 
concrete  in  a  pure  bending  mode  and  under  slow  load  conditions.  The  ultimate 
failure  strain  and  the  load  which  created  this  strain  are  then  used  to  ftstahlish 
paraiiKters  that  define  the  threshold  in  the  constitutive  model  that  will  be  used  for 
determining  if  cracking  will  occur  in  the  concrete  model.  Since  the  iiltiinaf<» 
strain  capacity  is  the  key  component  in  crack  prediction,  it  is  important  to 
accurately  determine  this  parameter. 

The  slow  load  beam  test  is  con^rised  of  a  simply  supported  beam  equally 
loaded  at  the  third  points.  Figure  22  shows  how  the  test  is  set  up  as  well  as  how 
the  dimensions  are  used.  The  beam  is  loaded  through  the  use  of  a  hydraulic  ram 
operated  under  electronic  control.  The  loading  of  the  beam  is  increased  on  a 
weekly  basis  until  the  beam  fractures. 


Rgure  22.  Geometry  and  loading  of  beam  used  in  slow  load  test 

Based  on  the  time  the  fracture  occurred  and  die  load  which  created  the 
fracture,  the  ultimate  stress  can  be  confuted  and  the  ultimate  strain  capacity  can 
be  measured.  The  threshold  of  cracking  used  in  the  constitutive  model  is  an 
interaction  of  strain  and  stress  as  shown  in  Figure  23.  The  ultimate  fracture  strain 
e^is  the  anchor  point  on  the  horizontal  axis  of  the  interaction  diagram  while  the 
ultimate  fracture  stress  oy  varies  with  time  as  a  function  of  the  modulus  of 
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Rgure  23.  Stress-strain  interaction  diagram  used  as  cracking  threshold  in 
constitutive  model  (Oj  =  stress  at  the  bottom  of  the  beam; 

E (r)  =  modulus  of  elasticity  at  time  r) 

elasticity.  In  the  figure  is  a  fixed  point  while  Oy^will  vary  with  time  since  it  is 
related  to  the  modulus  of  elasticity  E  through  the  following  equation: 

o,,  .  €,£«,)  (1) 

where  r,  is  the  time  elapsed  since  the  time  of  set  Therefore,  the  threshold  as 
presented  in  Figure  23  is  for  one  instant  in  time  and  will  be  different  for  various 
times.  For  example  will  be  less  for  times  earlier  in  the  life  of  the  concrete  and 
will  be  higher  for  times  later  in  the  life  of  the  concrete. 

The  ultimate  fracture  strain  cannot  be  obtained  directly  from  the  test  results 
since  it  is  the  strain  at  which  fracture  will  occur  when  there  is  no  stress  present. 
This  condition  will  not  likely  occur  in  actual  structures,  but  use  of  this  point  to 
anchor  the  curve  to  the  horizontal  axis  shown  in  Figure  23  is  necessary  to  identify 
the  curve  used  to  predict  failure  for  actual  conditions.  Such  a  case  is  the  slow 
load  beam  test  and  therefore  e^^must  be  extrapolated  from  the  results  obtained 
during  the  test.  When  the  fracture  occurs,  the  known  parameters  are  the 
magnitude  of  the  applied  load,  the  age  of  the  concrete,  and  the  strain  near  the 
bottom  of  the  beam  (taken  from  strain  gauges  mounted  on  the  beam).  Since  the 
center  of  the  beam  is  in  pure  bending,  the  stress  in  the  bottom  of  the  beam 
(modulus  of  rupture  can  be  computed  using  the  load  at  the  time  of  fracture 

with  the  following  equation: 
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ac 

II 


where 


Pf=  total  load  at  fracture 

a  =  distance  from  a  support  to  the  nearest  point  of  applied  load 
c  =  distance  from  the  bottom  of  the  beam  to  the  centroid 
I  =  moment  of  inertia  of  the  beam 


Since  the  time  of  fracture  is  known,  the  modulus  of  elasticity  at  the  time  of 
fracture  can  be  ascertained  from  the  test  data  obtained  during  matfirial  testing 

Prior  to  confuting  the  ultimate  fracture  strain  of  the  concrete,  the  slow  load 
fracture  strain  (strain  at  the  bottom  of  the  beam  €, ),  and  the  slow  load  fracture 
stress  (stress  at  the  bottom  of  the  beam  a, )  must  be  determined.  These  param¬ 
eters  are  not  the  value  of  stress  and  strain  at  actual  failure  but  at  90  percent  of 
failure  and  are  used  because  the  stress-strain  relationship  becomes  nonlinear  near 
failure  due  to  microcracldng  occurring  in  the  concrete.  Thrarefore,  the  slow  load 
fracture  stress  is  0.9  (p^.  A  linear  regression  analysis  is  performed  whidi  relates 
the  extreme  fiber  strain  to  the  ^plied  stress,  and  it  is  from  this  resulting  linear 
regression  that  the  slow  load  fracture  strain  can  be  determined  based  on  the  value 
ofo,. 

Once  o,,  E,  and  e,  have  been  determined  as  described  above,  the  nlriinatft 
fracture  strain  can  be  computed  by  using  the  following  equation: 


E(.0 


(3) 


where  t,  is  the  time  at  which  slow  load  fracture  stress  was  adiieved  in  the  test. 

Based  on  established  parametos  within  the  constitutive  model  and  the 
ABAQUS  model  parameters,  the  stress-strain  state  at  eadi  integration  point  in  the 
model  is  checked  against  the  cracking  criterion  (Figure  24)  to  determine  if  any 
(racking  is  predictecL  If  the  stress-strain  state  of  an  integration  point  exceeds  the 
(racking  criteion,  then  die  stifhiess  at  that  point  is  reduced  to  zero  in  the 
diretrion  of  the  maximum  prin(rpal  strain. 

The  stress-strain  interaction  diagram  is  also  used  to  provide  information  about 
integration  points  which  have  not  cracked.  The  constitutive  model  provides  a 
(racking  potential  of  uncracked  points  based  on  the  location  of  the  stress-strain 
state  on  the  diagram  relative  to  the  stress-strain  threshold  Output  for  (racking 
potential  is  provided  in  the  form  of  a  percaitage.  This  percentage  is  con:q)uted  by 
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dividing  the  distance  between  the  stress-strain  state  and  the  origin  (1 ,)  by  the 
distance  from  the  stress-strain  threshold  to  the  origin  (1  j  + 12)  (Figure  24).  The 
cracking  potential  is  valuable  information  since  it  will  indicate  whether  or  not  a 
structure  is  approaching  cracking  and  can  be  plotted  in  contours  to  indicate  areas 
where  there  may  be  a  high  potential  for  cracking  even  though  the  model  did  not 
actually  predict  cracking. 


Figure  24.  Cracking  potential  generation  for  a  specific  cracking  failure  surface 


Test 

The  test  setup  shown  in  Figure  22  is  based  on  ASTM  C  78-75  (ASTM  1992b). 
The  test  specimen  is  12  in.  wide.  The  loading  of  the  beam  was  based  on  applic¬ 
able  provisions  of  C  71-80  (USAGE  1949)  which  specify  a  loading  rate  that  will 
increase  the  stress  in  the  extreme  fiber  of  the  beam  by  25  psi  per  week.  Based  on 
the  dimensions  of  the  beam  used  for  the  testing  of  the  interior  mixture  for  the 
McAlpine  Lxx:k  Replacement  Project,  this  loading  rate  requires  that  a  load  of 
720  lb  be  applied  weekly  to  the  beam. 

Due  to  various  factors,  the  increase  in  the  loading  did  not  take  place  in  exactiy 
seven  day  increments  as  specified  in  ASTM  C  78-75  (ASTM  1992b).  Figure  25 
is  a  plot  of  the  load  versus  time  showing  how  the  prescribed  loading  varies  when 
compared  with  the  actual  loading.  The  actual  loading  will  be  used  for  compari¬ 
sons  in  this  study. 

The  beam  failed  at  a  total  load  PfOt  10,080  lb  with  concrete  109  days  old. 
Based  on  the  known  load  and  the  fact  that  in  the  center  of  the  beam  the  stress 
distribution  is  linear,  the  stress  in  the  bottom  of  the  beam  just  prior  to  failure  can 
be  computed  using  Equation  2.  The  distance  a  is  20  in.,  the  distance  to  the 
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Rgure  25.  Comparison  of  recommended  loading  of  slow  load  beam  and 
actual  applied  loading 


centroid  c  is  6  in.,  and  the  moment  of  inertia  /  is  1728  in'*.  Substituting  these 
values  results  in  the  following: 


o 


r 


(10,080  lb) (20  in.)  (6  in.)  . 

- — - -  =  350  psi 

2(1728  in.'*) 


(4) 


Since  the  modulus  of  rupture  is  350  psi,  the  slow  load  fracture  stress  will  be 
90  percent  of  this  value,  or  315  psi.  Based  on  a  315-psi  slow  load  fracture  stress, 
the  slow  load  fracture  strain  at  the  bottom  of  the  beam  e^can  be  obtained  frrom  the 
linear  regression  of  stress  and  strain  as  159.0  pm  and  is  reported  in  Bombich, 
Neeley,  and  Gamer.*  The  modulus  of  elasticity  of  the  concrete  can  be  conq)uted 
based  on  the  equation  used  in  the  constitutive  model.  Under  the  established 
loading,  the  age  of  the  concrete  for  a  stress  of  315  psi  would  be  97  days  and 
results  in  a  modulus  of  elasticity  at  day  97  of  4.18  x  10^  psi.  Based  on  the  above 
values,  e^can  be  computed  using  Equation  3  as  follows: 

6.  =  159.0  X  10'®  iniin.  +  - 315_psi -  ^  ^34 4  in./in.  (5) 

4.18  X  10®  psi 


'  Op.  cit 
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Since  the  input  for  ABAQUS  is  for  the  uniaxial  case,  the  strain  value  needed  is 
half  of  that  computed  or  1 17.2  x  lO"*  in./in.  and  represents  the  fracture  that  would 
be  obtained  from  a  direct  tensile  test 


Numerical  Analyses 

Numerical  analyses  were  performed  on  three  different  meshes.  Coarse, 
medium,  and  fine  meshes  were  used  and  are  shown  in  Figure  26.  The  coarse 
mesh  used  elements  that  were  4  in.  by  4  in.,  the  medium  mesh  used  elements  that 
were  3  in.  by  3  in.,  and  the  fine  mesh  used  elements  that  were  2  in.  by  2  in.  All 
three  meshes  modeled  only  the  60  in.  of  beam  located  between  the  supports,  so 
each  mesh  is  60  in.  by  12  in.  The  three  different  meshes  were  used  in  an  attempt 
to  determine  the  varying  effect  of  mesh  size  on  the  results. 
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Rgure  26.  Meshes  used  to  analyze  the  slow  load  beam  test 


Initially,  analyses  were  performed  on  each  of  the  three  meshes  using  half  of 
the  ultimate  firacture  strain  e^-of  1 17.2  in./in.  This  resulted  in  a  predicted  failure 
at  three  different  times  of  the  analysis  as  shown  in  Table  3.  As  should  be 
expected,  the  coarse  mesh  predicted  the  latest  day  of  failure  since  all  three 
meshes  assumed  the  same  ultimate  fracture  strain.  Two  factors  affect  this 
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Table  3 

Predicted  Time  of  Failure  of  Slow  Load  Beam  Test  Based  on 
h/2e,=  117.2  in./in.  


Mesh  Density 

Coarse 

Medium 

Rne  B 

Predicted  day  of  failure 

90.5 

86.0 

81.5  1 

behavior.  The  distribution  of  strain  is  linear  from  the  top  to  the  bottom  of  the 
beam,  and  integration  points  are  farther  away  from  the  extreme  tensile  fiber  for 
coarse  meshes  when  compared  with  fine  meshes.  Therefore,  it  will  require  more 
load  for  the  coarse  mesh  to  attain  the  level  of  strain  that  occurs  in  the  medium  and 
fine  meshes  and  likewise  when  comparing  the  medium  mesh  to  the  fine  mesh. 
Regardless  of  the  mesh  chosen,  the  resulting  predicted  failure  occurs  earlier  than 
the  actual  failure  in  the  laboratory  in  all  three  cases,  which  indicates  a 
conservative  application  of  the  failure  criteria. 

In  order  to  assess  the  failures  predicted  by  each  of  the  three  nieshes,  the 
analyses  were  performed  again  for  each  mesh  except  that  the  value  used  for  half 
of  the  ultimate  fracture  strain  was  adjusted.  This  adjustment  conqwnsated  for  the 
location  of  the  integration  point  with  respect  to  the  bottom  of  the  beam  so  that 
when  the  strain  at  the  extreme  fiber  reached  the  actual  one-half  of  value,  the 
strain  at  the  integration  point  location  would  correspond  to  the  fracture  strain. 

The  adjusted  values  for  one  half  of  were  computed  based  on  the  assunq)tion  of 
linear  strain  distribution  and  by  using  similar  triangles.  Table  4  shows  the 
adjusted  values  of  one  half  Cy^used  in  the  analyses  and  the  predicted  day  of  failure 
for  each  mesh.  As  might  be  expected,  all  three  meshes  predict  failure  of  the  beam 
on  the  same  day.  As  in  the  previous  case  though,  the  failure  occurs  at  much 
earlier  times  in  the  analyses  than  in  the  laboratory. 


Table  4 

Predicted  Time  of  Failure  of  Slow  Load  Beam  Test  Based  on 
Adjusted  Values  of  1/26,  in. 


Mesh  Density 

Coarse 

Medium 

Rne  1 

Adjusted  1/2e,(ln7in.) 

100.7 

104.8 

108.9  1 

Predicted  day  of  failure 

74.5 

74.5 

74.5  I 

Time-histories  of  strain  measurements  taken  from  two  locations  on  the  beam 
during  the  test  are  plotted  in  Figures  27  and  28.  Also  plotted  in  Rgure  27  are  (he 
predicted  total  strain  histories  from  the  numerical  analyses  which  used  a  constant 
value  of  117.2  in./in.  for  l/2€y  for  all  three  meshes,  while  Figure  28  includes 
predicted  strain  histories  for  the  analyses  which  used  the  adjusted  values  of  l/2ey 
as  shown  in  Table  4.  As  can  be  seen,  the  predicted  strains  from  the  analyses 
correlate  very  well  with  the  actual  strains  measured.  The  strain  predicted  by  the 
fine  mesh  matches  closely  the  strain  that  was  measured  at  the  bottom  of  the  beam. 
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Ftgure  27.  Strain  time-histories  of  measured  strains  in  actual  test  and  predicted 
strains  from  numerical  analyses  using  constant  value  of  1 17.2  inTin. 
fore, 

while  the  results  from  the  medium  and  coarse  meshes  fall  in  between  the  two  sets 
of  measured  strains.  This  is  due  to  the  location  of  the  integration  points  in  each 
of  the  analyses.  The  bottom  integration  points  of  the  fine  mesh  are  only  0.423  in. 
from  the  base  of  the  beam,  while  the  bottom  integration  points  for  the  n^dium 
mesh  and  the  coarse  mesh  are  0.634  in.  and  0.845  in.  from  the  bottom  of  the 
beam,  respectively.  This  puts  the  location  of  the  integration  points  for  the 
medium  and  coarse  meshes  close  to  die  middle  point  between  the  two  actual 
gauges.  Figures  27  and  28  also  demonstrate  the  effects  of  using  a  common  value 
for  l/2€^as  opposed  to  using  an  adjusted  value.  When  a  single  value  is  used,  the 
analyses  (Figure  27)  all  fail  at  different  times  but  the  strain  value  at  which  they 
fail  is  approximately  the  same.  However,  when  the  values  are  adjusted  for  the 
mesh  (Figure  28),  the  failure  for  each  analysis  occurs  on  the  same  day  but  the 
strain  values  at  which  failure  occurs  vary. 

Questions  may  be  raised  concerning  the  fact  that  the  predicted  strains  match 
the  actual  measured  strains  so  closely,  yet  the  failure  times  are  so  different.  The 
discrepancy  occurs  because  the  strain  used  by  the  constitutive  model  to  predict 
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Rgure  28.  Strain  time-histories  of  measured  strains  in  actual  test  and  predicted 
strains  from  numerical  analyses  using  adjusted  values  for 

fracture  is  different  than  the  total  strain  that  is  plotted  in  Figures  27  and  28.  The 
strain  plotted  in  these  figures  includes  shrinkage  strains  due  to  the  autogenous 
shrinkage  which  is  included  in  the  constitutive  model.  These  strains  are  not  part 
of  the  fracture  strain  used  to  predict  cracking  since  these  strains  do  not  produce  a 
mechanical  type  of  strain  which  causes  fracture.  The  shrinkage  strains  can  be 
compared  with  thermal  strains  which  can  create  volumetric  changes  without 
putting  a  body  under  stress  if  there  are  no  restraining  conditions  present 
Therefore,  shrinkage  strains  associated  with  autogenous  shrinkage  wiU  not 
contribute  to  the  cracking  of  the  beam  since,  for  a  simple  beam  arrangement, 
there  is  nothing  restraining  the  movement  associated  with  these  strains 
Comparing  Figure  28  with  Figure  29  shows  how  the  total  strain  is  different  from 
the  fracture  strain.  As  can  be  seen  in  Figure  29,  the  strains  are  higher  than  those 
in  Figure  28  since  the  shrinkage  strains  have  been  removed. 

As  can  be  seen  in  Figure  29,  the  fracture  strains  at  failure  are  approximately 
the  same  values  as  the  maximum  measured  strain  before  failure  at  the  extreme 
fiber  (second  curve  in  Figure  29).  This  indicates  that  the  failure  strain  in  the 
constitutive  model  accurately  reflects  the  strain  state  at  failure.  However,  while 
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Rgure  29.  Strain  histories  of  measured  strains  in  actual  test  and  predicted 

cracking  strains  from  numerical  analyses  using  adjusted  values  of  e, 
(similar  to  Rgure  28)  (N.A.  =  neutral  axis) 

the  strain  states  may  be  similar,  the  stress  states  are  different.  The  load  at  which 
the  beam  failed  in  the  numerical  studies  is  approximately  30  percent  less  than  the 
load  at  which  the  actual  beam  failed.  This  indicates  that  the  strain  to  which  the 
model  is  calibrated  should  be  adjusted  to  remove  shrinkage  strains  prior  to 
performing  the  calibration  in  order  to  capture  the  correct  behavior. 

During  the  study,  another  factor  came  into  consideration  regarding  the  creep. 
Currently  the  constitutive  model  includes  a  temperature  effect  for  the  creep.  This 
temperature  effect  was  established  for  high  ten^rature  applications  such  as 
concrete  used  in  nuclear  reactors  and  for  which  there  were  data  available. 
However,  the  same  relationships  are  used  in  the  constitutive  model  used  in  NISA 
studies,  and  physical  data  do  not  exist  to  support  the  use  of  such  a  temperature 
effect  in  the  constitutive  model  for  concrete  in  the  temperature  ranges  typically 
observed  in  mass  concrete  structures.  This  temperature  effect  can  be  seen  when 
the  predicted  strain  histories  shown  in  Figure  30  are  compared  with  those  in 
Figure  28.  The  strains  shown  in  Figure  30  are  a  result  of  performing  analyses  in 
which  the  initial  temperature  is  not  specified  in  the  ABAQUS  input  deck,  and 
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Rgure  30.  Strain  time-histories  of  measured  strains  in  actual  test  and  predicted 
strains  from  numerical  analyses  using  adjusted  values  for  e,and  no 
initial  temperature  of  70  °F 

therefore  a  temperature  of  0  °F  will  be  used  by  ABAQUS  to  perform  the  analysis. 
This  0  °F  assumption  will  be  transferred  to  the  constitutive  model  and  the 
resulting  change  is  seen  in  Figure  30.  The  original  analyses  used  an  initial 
temperature  of  70  °F,  for  which  results  are  presented  in  Figure  28,  which 
corresponds  to  the  temperature  used  to  perform  the  test 

As  can  be  seen  in  Figure  30,  the  predicted  strains  for  the  case  where  the 
tenqwrature  is  not  initialized  at  70  °F  are  much  lower  than  those  predicted  when 
the  70  °F  initial  temperature  is  used.  Analyses  of  all  three  meshes  produced 
failure  at  day  109  of  the  analysis  which  corresponds  to  the  actual  beam  failure  in 
the  laboratory.  While  this  would  seem  to  indicate  that  the  case  using  no  initial 
temperature  accurately  captures  the  beam  behavior,  review  of  the  predicted  strain 
histories  shows  how  this  case  underpredicts  the  actual  strain.  That  this  case 
matches  the  same  day  of  failure  as  the  actual  tests  is  at  least  partly  due  to  the  fact 
that  the  fracture  strain  does  not  include  shrinkage  strains.  Plotting  the  fracture 
strains  for  this  case  shows  that  these  strains  follow  more  closely  the  measured 
strains.  This  is  demonstrated  in  Figure  31  where  the  predicted  strains  (the 
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Rgure  31 .  Strain  time-histories  of  measured  strains  in  actual  test  and  predicted 
strains  from  numerical  analyses  using  adjusted  values  for  e,and  no 
initial  temperature  of  70  °F  but  with  shrinkage  strains  included 
(similar  to  Figure  30) 


fracture  strains)  closely  approximated  the  measured  strains  located  4.5  in.  from 
the  neutral  axis. 


Conclusion 

The  predictions  made  by  the  numerical  analyses  of  the  slow  load  beam  test 
showed  good  agreement  with  the  measured  data  taken  during  the  test  While  the 
comparison  of  total  strains  was  good,  the  numerical  studies  showed  that  the  beam 
failed  much  earlier  in  the  loading  history  than  actually  occurred.  This  indicates 
that  the  method  of  calibration  must  be  adjusted.  Since  it  is  true  that  strains  in  a 
single  beam  due  to  autogenous  shrinkage  will  not  induce  cracking,  these  strains 
should  not  be  included  in  the  fracture  strain  used  to  predict  cracking.  However, 
the  strain  that  is  being  measured  in  the  beam  test  is  the  total  strain  and  does 
include  all  strains  that  are  occurring  in  the  beam,  including  any  shrinkage  strain. 


36 


Chapter  4  Slow  Load  Beam  Test 


Therefore,  while  the  shrinkage  strain  does  not  contribute  to  the  cracking  of  the 
beam,  it  is  still  present  and  must  be  taken  into  account  when  comparing  with 
strains  measured. 

Strains  associated  with  creep  should  be  considered  in  a  similar  manner  as  the 
shrinkage  strains.  Since  creep  is  the  elastic  deformation  of  the  concrete,  it 
generally  will  not  contribute  to  cracking  of  the  concrete  imless  some  sort  of 
restraint  is  introduced.  Therefore,  the  correct  method  for  modeling  the  strain 
would  be  to  separate  the  creep  strain  from  the  total  strain  used  to  evaluate  fracture 
as  is  done  with  the  strain  associated  with  autogenous  shrinkage.  The  proprietary 
owner,  ANATECH,  of  the  ANACUP-U  software  (1992b)  which  contains  the 
constitutive  model,  is  in  the  process  of  updating  the  model  to  separate  the  creep 
strain  from  the  total  strain  which  should  correct  the  discrepancies  observed  in  the 
above  study.  Once  this  update  has  occurred,  accounting  for  the  creep  strain  prior 
to  calibration  will  also  need  to  be  considered.  In  addition,  ANATECH  currently 
has  plans  to  remove  the  temperature  dependency  of  the  creep  relationship  on 
tenq)erature  since  there  are  no  data  available  to  support  this  dependency  at 
temperatures  of  100  °F  and  less. 

Another  matter  which  should  be  considered  in  future  studies  is  a  possible 
reduction  of  the  failure  strain  input  into  the  ABAQUS  deck.  The  computation  of 
the  failure  strain  as  currently  performed  uses  stress  and  strain  occurring  at  the 
extreme  tensde  fiber  of  the  beam.  It  was  demonstrated  that  the  coarser  the  mesh, 
(he  longer  it  takes  to  actually  attain  fracture.  While  this  should  be  intuitively 
obvious,  it  does  point  out  that  consideration  should  be  given  to  the  fact  that  the 
numerical  analysis  does  not  perform  its  calculation  at  the  extreme  fiber  and  that 
the  cracking  threshold  should  be  adjusted  to  account  for  this  discrepancy. 
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Introduction 

The  rapid  load  beam  test  is  used  to  measure  the  fracture  strain  in  a  rapid 
loading  environment  to  determine  strains  at  fracture  without  the  inclusion  of 
time-dependent  factors  such  as  creep  and  shrinkage.  The  results  of  the  rapid  load 
beam  test  are  not  required  for  input  into  the  constitutive  model,  but  it  will  be  of 
value  to  make  conq)arisons  between  rapid  load  beam  test  data  and  results 
predicted  in  a  numerical  analysis  utilizing  the  constitutive  model. 

The  setup  of  the  rapid  load  beam  test  is  the  same  as  the  slow  load  beam  test; 
i.e.,  it  is  comprised  of  a  sinq)ly  supported  beam  equally  loaded  at  the  third  points. 
Figure  22  shows  how  the  test  is  set  up  as  well  as  the  dimensions  used.  The  beam 
is  loaded  through  the  use  of  a  hydraulic  ram  operated  under  electronic  control. 
The  loading  of  the  beam  is  increased  continuously  at  a  rapid  rate  until  the  beam 
fractures. 


Test 

The  test  setup  shown  in  Figme  22  is  based  on  ASTM  C  78-75  (ASTM  1992b). 
The  test  specimen  is  12  in.  wide.  The  loading  of  the  beam  was  based  on 
applicable  provisions  of  C  71-80  (USAGE  1949)  which  specify  a  loading  rate  that 
will  increase  the  stress  in  the  extreme  fiber  of  the  beam  by  40  psi  per  minute. 
Based  on  the  dimensions  of  the  beam  used  for  the  testing  of  the  interior  mixture 
for  the  McAlpine  Lock  Replacement  Project,  this  loading  rate  requires  that  a  load 
of  1,152  lb  be  applied  each  minute  to  the  beam.  Unlike  the  slow  load  beam  test, 
the  rapid  load  beam  test  is  set  up  for  automatic  loading.  The  loading  takes  place 
as  specified  and  readings  of  load  and  strain  are  recorded  every  ten  seconds. 

Two  rapid  load  beam  tests  were  performed;  one  when  the  concrete  was  7  days 
old,  and  the  other  at  day  109  when  the  slow  load  beam  failed.  Tlie  test  performed 
at  day  7  failed  with  a  load  of  5,509  lb  which  meant  die  stress  in  the  extreme 
tensile  fiber  was  191.3  psi.  For  the  beam  tested  at  day  109,  the  beam  failed  under 
a  load  of  15,414  lbs  which  induced  a  stress  of  535.2  psi.  Maximum  tensile 
strains  at  failure  were  63  ^m  for  the  day  7  test  and  1 19  fj-m  at  day  109. 
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Numerical  Analyses 

As  in  the  slow  load  beam  test,  numaical  analyses  were  performed  on  three 
different  meshes.  A  coarse,  medium,  and  fine  mesh  were  used  and  are  shown  in 
Figure  26.  The  coarse  mesh  used  elements  that  were  4  in.  by  4  in.,  the  mpyiinm 
mesh  used  elements  that  were  3  in.  by  3  in.,  and  the  fine  mesh  used  elp-niftnts  that 
were  2  in.  by  2  m.  Analyses  of  the  rapid  load  test  were  performed  at  7  days  and 
at  109  days  to  match  the  tests  performed  in  the  laboratory.  For  both  analyses,  the 
time  history  of  the  beam  included  the  time  up  until  loading.  The  beam  was 
modeled  as  being  fully  supported  along  the  bottom  of  the  beam  up  until  the  Hmp. 
of  loading  for  both  analyses. 

Results  of  the  analyses  conqjared  with  the  actual  test  data  are  shown  in 
Figures  32  and  33  in  the  form  of  stress-strain  plots.  The  numaical  results  shown 
have  all  been  adjusted  so  that  the  stresses  and  strains  at  the  beginning  of  loading 
are  zero.  In  addition,  the  resulting  strains  reported  from  the  analyses  at  the 
integration  points  have  been  interpolated  so  that  strains  at  the  extreme  fibors  are 
actually  plotted.  This  was  done  since  the  strains  obtained  from  the  test  Hata  had 
also  been  adjusted  to  the  extreme  fibers.  As  a  result  of  these  adjustments,  the 
curves  of  the  three  sets  of  numerical  data  all  fall  on  top  of  one  another  and  give 
the  appearance  of  only  one  set  of  munerical  data. 

Figure  32  shows  the  results  of  the  rapid  load  beam  test  at  an  age  of  7  days 
compared  with  those  of  the  numerical  analyses  performed.  While  there  is  some 
discrq)ancy  between  the  two  sets  of  data,  the  difference  is  not  large.  Since  the 
noodulus  of  elasticity  can  be  determined  from  stress-strain  plots,  some  explana¬ 
tion  of  the  difference  may  be  discovered  from  conqruting  the  predicted  modulus 
based  on  the  curves.  Based  on  the  strain  values  at  an  applied  stress  of  17.3  psi 
and  1 18.7  psi,  the  modulus  for  the  actual  test  is  approximately  3,343,000  psi, 
while  using  strains  associated  with  stresses  of  20  psi  and  240  psi  for  the  numeri¬ 
cal  studies  results  in  a  modulus  of  approximately  2,37 1,000  psi.  The  modulus 
predicted  by  the  numerical  study  matches  well  with  the  modulus  of  elasticity 
curve  used  in  the  model  v(iiidi  has  a  7-day  modulus  of 2,404,000  psi.  The 
numerical  results  also  compare  favorably  with  that  reported  in  Bombich,  Neeley, 
and  Gamer'  whidi  shows  a  modulus  of  2,300,000  psi  at  day  7.  From  this  it  is 
evident  that  the  numerical  model  is  predicting  the  stress-strain  response  correctly 
based  on  the  information  for  which  it  was  calibrated.  While  measuring  errors  or 
changes  in  mixture  proportions  are  possible  reasons  for  the  differences  seen,  they 
are  not  likely  to  occur.  A  more  plausible  reason  for  the  disaepancy  between  the 
results  predicted  by  the  numaical  study  and  the  results  from  the  tests  can  be 
attributed  to  the  fact  that  the  concrete  used  in  the  test  to  determine  the  modulus  of 
elasticity  is  wet-sieved  prior  to  testing  while  the  conaete  for  the  beam  test  is  not 
The  fact  that  the  concrete  used  in  the  beam  test  is  not  wet-sieved  means  that  it 
contains  more  aggregate.  Smce  the  aggregate  has  a  mudi  higher  mnduhis  than 
the  concrete  at  early  ages  (i.e.,  7  days),  the  presence  of  more  aggregate  in  a 
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Figure  32.  Measured  strain  versus  predicted  strain  for  rapid  load  test  at 
concrete  age  of  7  days 


conaete  mixture  will  result  in  a  stiffer  conaete  due  to  bridging  that  will  occur 
between  aggregate  particles. 

It  should  be  noted  that  while  the  numerical  model  appears  to  have  performed 
well  in  capturing  the  stress-strain  behavior  of  the  conaete,  it  spears  that  at  early 
rimfis  it  may  overpredia  the  aacking  capacity  of  the  conaete.  The  numerical 
study  indipatpri  the  strain  in  the  conaete  was  nearly  double  of  that  tested,  and  the 
stress  was  50  percent  greata.  E>ue  to  the  diffaence  in  the  modulus  of  the  actual 
beam  test  and  the  nnodulus  used  in  the  constitutive  model,  determining  the  cause 
for  these  diffaences  is  difficult  The  diffaence  could  be  a  result  of  the  constitu¬ 
tive  model  being  unable  to  propaly  capture  the  behavior  of  the  conaete  but  this 
cannot  be  expressed  with  certainty  unless  the  constitutive  model  is  calibrated  to 
the  conaete  mixture  actually  being  used  in  the  test  Thaefore,  furtha  tests  are 
needed  to  eitha  confirm  or  invahdate  the  observations  made  in  this  study. 

Figure  33  shows  the  conq)arison  between  the  rapid  load  test  performed  at  day 
109  and  the  numerical  studies  performed.  Conq)arison  of  test  results  with 
numerical  results  are  very  good  in  this  plot  Afta  the  initial  loading,  the  two 
curves  look  to  be  almost  completely  parallel.  This  is  evidenced  in  the  calculation 
of  the  modulus  of  each  curve.  For  the  actual  test  data  using  strains  with 
associated  stresses  of  100  psi  and  442  psi,  a  modulus  of  4,343,000  is  calculated, 
while  the  numerical  studies  strains  with  associated  stresses  of  40  psi  and  440  psi 
result  in  a  modulus  of  4,326,000  psi.  The  cahbration  of  the  numerical  model 
shows  that  at  day  109  the  modulus  of  the  material  is  4,380,000  psi  which 
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Figure  33.  Measured  strain  versus  predicted  strain  for  rapid  toad  test  at 
concrete  age  of  109  days 


compares  well  with  the  modulus  values  coiiq)uted  based  on  test  data  and 
numerical  data.  There  is  no  actual  data  for  the  modulus  of  elasticity  in  the 
materials  rqKut^  but  if  an  intopolation  between  the  values  givai  at  day  90  and 
day  180  is  performed  the  resulting  modulus  is  4,180,000  psi,  whidi  once  again 
con^ares  favorably  with  both  sets  of  data.  The  similarity  of  the  modulus  test  and 
beam  test  at  this  point  in  time  can  be  attributed  to  the  fact  that  the  values  of  the 
rock  noodulus  and  concrete  modulus  are  much  closer  at  day  109  than  at  day  7. 

Unlike  the  rapid  load  test  performed  at  day  7,  the  stress  and  strain  values  at 
failure  conq)are  very  well.  The  strain  at  fracture  is  nearly  identical  (119  jum),  and 
the  applied  stress  is  only  slightly  higher  for  the  experimoital  results  compared 
with  the  nummcal  results.  These  results  indicate  that  the  material  model  at  later 
times  predicts  the  behavior  of  the  concrete  very  well,  including  the  load  and 
strain  at  which  it  fails. 


Conclusions 

The  results  of  the  num^cal  studies  of  the  r<^id  load  beam  test  conq)ared 
favorably  with  the  results  of  the  actual  experimental  results.  The  differences  in 
the  experimental  results  and  numerical  results  were  more  than  desired,  but  it  is 
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unclear  if  this  is  an  actual  (Ufference  or  not  since  the  modulus  of  elasticity 
predicted  in  the  experimental  results  does  not  match  the  expected  modulus  for  the 
conCTete  at  an  age  of  7  days.  With  respect  to  the  design,  the  model  performed 
extremely  well;  i.e.,  the  modulus  of  elasticity  values  predicted  by  the  stress-strain 
curves  matched  the  values  to  which  the  model  was  calibrated.  Since  questions 
still  remain  about  the  strains  at  early  times,  further  testing  should  be  performed 
which  will  permit  some  of  the  questions  raised  concerning  early  time  strains 
during  rapid  loading  to  be  answered. 
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6  Conclusions  and 
Recommendations 


The  data  contained  within  this  report  do  not  answCT  all  of  the  questions 
regarding  the  validity  of  the  constitutive  model  used  in  a  NISA,  but  they  do 
provide  insight  into  where  the  model  does  and  does  not  paform  weU.  Overall, 
the  results  presmted  indicate  that  the  constitutive  model  is  capable  of 
^proximating  results  of  actual  tests.  Presented  below  are  specific  conclusions 
which  can  be  made  regarding  the  constitutive  model  as  weU  as  recommendations 
of  vdiat  may  be  needed  in  the  future  to  assist  in  validation  of  the  model. 


Conclusions 

The  analyses  performed  for  conqtarison  with  multiaxial  creep  data  conqtared 
very  well  with  test  data  for  low  levels  of  compressive  loadings  but  did  not 
paform  well  for  high  con^ressive  loads  outside  of  the  elastic  creep  range  for 
■wiiich  the  model  was  intended.  The  fact  that  the  model  cannot  c^ture  the 
behavior  of  aeep  imder  high  con^tressive  loads  does  not  detract  firom  the  ability 
of  the  model  to  be  used  in  analyzing  massive  conorete  structures  since 
conq)ressive  loads  will  rarely  exceed  loadings  greater  than  40  percent  of  the 
con5>ressive  strength,  whidi  is  the  basis  of  the  model.  It  is  possible  that  this  limit 
might  be  exceeded  in  dynamic  load  situations,  so  caution  should  be  exercised  if 
attempting  to  use  the  model  for  this  type  of  condition.  Based  on  the  results  of  the 
studies  reported  in  Chapter  2  for  loadings  under  40  percent  of  the  con5)ressive 
load,  the  creep  curves  which  are  based  on  uniaxial  tests  do  provide  reasonable 
predictions  for  other  states  of  stress. 

The  analyses  for  the  creep  strain  recovery  were  not  as  successfiil  as  the 
multiaxial  creep  studies.  Results  of  the  analyses  did  not  compare  well  with  the 
test  data.  As  noted  in  Chapte  3,  the  constitutive  model  was  developed  under  the 
assunqjtion  of  a  linear  creep  curve  which  allows  the  principle  of  superposition  to 
be  used.  However,  the  principle  of  superposition  does  not  accurately  predict  the 
deCTease  in  strains  due  to  unloading  (Bazant  1988).  This  would  explain  the  large 
differences  between  experimental  and  nummcal  results  observed  in  Chapter  3. 
The  cases  observed  are  for  instances  where  the  imloading  occurs  suddenly  and 
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the  CTeep  recovery  is  overpredicted.  For  cases  in  which  the  unloading  occurs 
gradually,  such  as  in  massive  concrete  structures,  the  model  may  be  adequate. 

Strain  con5)arisons  of  total  strain  were  very  good  for  the  slow  load  beam  test 
The  strain  histories  predicted  by  the  three  different  meshes  used  in  the  numerical 
studies  fell  between  the  strain  histories  recorded  by  the  gauges  located  at  the 
bottom  of  the  beam  and  at  1-1/2  in.  from  the  bottom  of  the  beam.  This  indicates 
that  the  numerical  predictions  were  very  good  since  the  integration  points  from 
which  the  numerical  data  were  taken  are  located  between  the  points  where  the 
strain  gauges  on  the  beam  were  mounted.  While  the  total  strains  compared  well, 
the  load  at  w^ch  fracture  occurred  was  significantly  different.  The  numerical 
studies  predicted  a  failure  of  the  beam  at  tpproximately  70  percent  of  the  load 
vi^ch  caused  failure  in  the  laboratory.  This  difference  is  attributed  to  the  fact 
that  the  constitutive  model  subtracts  strain  associated  with  autogenous  shrinkage 
from  total  strain.  It  is  true  that  strains  associated  with  autogenous  shrinkage  do 
not  contribute  to  fracture;  however,  the  model  is  currently  calibrated  to  the  total 
strain  as  measured  during  the  test.  This  means  that  the  current  method  of 
calibration  creates  a  situaticMi  where  a  model  that  excludes  shrinkage  strain  as 
part  of  the  fracture  strain  is  calibrated  to  a  strain  that  does  include  shrinkage 
strain.  This  is  of  course  inconsistent,  and  the  strains  obtained  during  testing 
should  be  adjusted  so  that  the  constitutive  model  can  be  properly  calibrated. 

Comparisons  made  for  the  npid  load  beam  test  differed  for  the  two  times  at 
which  it  was  performed.  The  test  performed  at  an  age  of  7  days  resulted  in  a 
difference  in  the  total  strains  reported,  but  due  to  the  fact  that  the  stress-strain 
curve  plotted  based  on  the  numerical  studies  provided  a  modulus  of  elasticity 
value  close  to  that  expected  for  the  conaete  tested  at  an  age  of  7  days,  it  can  be 
reasoned  that  wet-sieving  the  concrete  has  a  definite  effect  on  the  modulus  of 
elasticity  at  early  times.  At  day  109,  \vfren  the  other  rapid  load  beam  test  was 
paformed,  agreement  betweoi  the  measured  results  and  numerical  results  was 
very  good  and  indicated  that  the  model  can  accurately  predict  strains  in  a  rapid 
load  situation. 


Recommendations 

Based  on  studies  perfonned  on  multiaxial  creep  tests,  it  appears  that  the 
calibration  of  the  creep  function  as  it  is  currently  performed  is  sufficient  for  use 
in  NISA  studies  as  they  are  currently  conducted.  The  only  recommendation 
regarding  the  aeep  model  is  that  caution  should  be  exercised  in  using  the  current 
mcxlel  in  a  dynamic  type  analysis  since  it  will  not  paform  well  for  cuses  of 
nonlinear  creq)  and  for  stress  states  that  are  high  compared  with  the  compressive 
strength  of  the  concrete. 

For  the  studies  on  creep  recovery  it  is  clear  that  a  linear  definition  of  creep 
docs  not  sufficiently  cr^ture  the  strain  recovery  that  occurs  when  sudden  drops  in 
load  exxur.  The  data  available  do  not  allow  an  assessment  to  be  made  regarding 
creep  strain  recovery  when  the  load  is  reduced  gradually.  It  is  possible  that  the 
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model  could  properly  predia  creep  strain  recovery  und^  gradual  unloadings  as 
occur  in  a  massive  concrete  structure,  but  confirmation  of  this  is  needed.  It  is 
therefore  recommended  that  at  least  one  series  of  creep  tests  be  performed  which 
allow  a  gradual  unloading  of  the  specimens  and  another  series  which  would 
introduce  the  unloading  in  a  series  of  sudden  drops  in  load.  Tests  such  as  these 
would  provide  the  information  necessary  for  determining  the  accuracy  of  the 
constitutive  model  with  respect  to  c^turing  creep  recovery  strain 

The  conqtarisons  made  of  total  strains  from  the  numerical  studies  provided 
excellent  agreement  with  the  laboratory  results.  However,  evaluations  of  the 
strains  made  it  r^paroat  that  the  calibration  was  not  being  performed  properly. 
Currently,  the  fracture  strain  is  calibrated  to  the  measured  strain  from  the 
laboratory  tests,  whicih  is  the  total  strain.  In  its  cxnrent  ftjrm,  the  fracture  strain  is 
the  total  strain  reduced  by  the  strain  asscxaated  with  autogenous  shrinVagp.  it  is 
clear  that  for  the  fracture  strain  to  be  properly  calibrated,  the  strains  measured  in 
the  laboratory  must  be  reduced  by  an  estimated  amount  of  shrinlcagp  strain  prior 
to  pofoiming  the  c:alibration.  As  noted  above,  a  future  version  of  the  constitutive 
model  will  remove  the  strains  associated  with  cxeq)  as  well.  It  is  recommended 
that  in  the  future,  prior  to  calibrating  the  fracture  strain  as  outlined  in  ETL  1 1 10- 
2-365  (USAGE  1994),  the  strain  at  fracture  obtained  from  the  slow  load  beam  test 
be  reduced  by  the  amount  of  strain  predicted  by  the  calibrated  shrinlcagp  curve  in 
the  constitutive  mcxlel.  This  will  ensure  that  the  cx>nq)arison  is  of  frac^ture  strain 
rather  than  total  strain. 

While  the  r^id  load  beam  test  conqiarisons  were  very  good  at  a  concrete  age 
of  109  days,  the  disaqiancy  between  numerical  results  and  test  results  at  the 
conaete  age  of  7  days  was  inconclusive  because  the  modulus  of  elasticity 
calculated  from  the  stress-strain  curve  for  the  test  results  did  not  matrh  the 
expected  value.  It  is  recommended  that  an  additional  beam  test  be  performed  at 
an  age  of  7  days  and  that  tests  for  measuring  the  modulus  of  elasticity  be 
paformed  on  the  same  concrete  mixture  used  in  the  beam  tests.  This  will  likely 
require  that  largo-  cyiinders  than  are  normally  used  due  to  the  presence  of  larger 
size  aggregate  because  no  wet-sieving  will  be  performed.  It  is  also  recommended 
that  these  tests  be  performed  using  the  interior  ccmcrete  mixture  from  the 
McAlpine  L(x;k  Replacement  Project.  This  will  allow  confirmation  of  the  ability 
of  the  constitutive  model  to  accurately  define  the  behavior  due  to  r^id  loads  at 
young  ages  or  to  identify  an  area  that  may  need  to  be  revised  in  the  constitutive 
model. 
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